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Abstract 

After a brief presentation of the MACRO detector we discuss the data on atmospheric neutrinos 
and neutrino oscillations, on high energy (E v > 1 GeV) neutrino astronomy, on indirect searches for 
WIMPs and low energy (E v ^7 MeV) stellar collapse neutrinos. 



i/-) 1 Introduction 

in 

MACRO is a multipurpose underground detector designed to search for rare events in the cosmic radiation. 
It was optimized to search for the supermassive magnetic monopoles (MMs) predicted by Grand Unified 
Theories (GUT) of the electroweak and strong interactions. MACRO performs measurements in areas 
of astrophysics, nuclear, particle and cosmic ray physics. In this paper we shall discuss atmospheric 
neutrinos and neutrino oscillations, high energy (E v ^ 1 GeV) neutrino astronomy, indirect searches for 
WIMPs and searches for low energy £ 7 MeV) stellar collapse neutrinos. The detector has global 
dimensions of 12 x 9.3 x 76.5 m 3 and provides a total acceptance to an isotropic flux of particles of 
^H. ~ 10, 000 m 2 sr. The total mass is ~ 5300 t. A cross section of the detector is shown in Fig. [j] jp. It 

has three sub-dctectors: liquid scintillation counters, limited streamer tubes and nuclear track detectors. 
The mean rock depth of the overburden is ~ 3700 m.w.e.; the minimum is 3150 m.w.e. This defines 
the minimum downgoing muon energy at the surface as ~ 1.3 TeV in order to reach the underground 
lab. The average residual energy and the muon flux at the lab depth are ~ 310 GeV and ~ 1 m~ 2 h _1 , 
H ' respectively. 

2 Atmospheric v. Neutrino oscillations. 

Neutrino oscillations. For massive neutrinos one has to consider the weak flavour eigenstates (v e , 
v T , which are relevant for 7r — > fiv^, K— ► uv^ and /i — > ev e v^ decays and for v interactions) and mass 
eigenstates (vx, v-i, ^3, relevant for propagation). The flavour eigenstates are linear combinations of the 
mass eigenstates, vi = Y] Ui m v m . For two flavours 

= V2COS6 + v^sinO . . 

i T = —visinB + v^cosO 

where 9 is the mixing angle. The rate of disappearance is given by 

P(v M -> Ufj) ~ 1 - sin 2 29 sin 2 '(1.27 Am 2 L/E u ) (2) 

Am 2 = 77j 2 „ 3 — m v 2 , L = distance (in m) from the decay point to the interaction point. 

Only disappearance experiments have been performed until now; the final proof of ^-oscillations will 
come from appearance experiments, P(f^ —> v T ) = 1 — P(v^ — > 1^), in which one observes neutrinos 
which have not been produced in the 7T, K decays. 
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Figure 1: Event topologies induced by i/ M interactions in and around MACRO. The black boxes are scintillator 
hits. The track direction and versus are measured by the streamer tubes and time of flight (ToF). 



Atmospheric neutrinos. Primary cosmic rays produce in the upper atmosphere pions and kaons, which 

decay, ir — > [iv, K — » /Ltf; and /i* — ► v ^ + v e +e ± . Neutrinos are produced in a spherical surface, ~ 
10 km from the Earth, and move towards the Earth. 

Neutrino-induced Upward Going Muons. They are identified using the streamer tube system (for 
tracking) and the scintillator system (for time-of- flight (ToF) measurement). A rejection factor of at least 
10 5 is needed in order to separate the up-going muons from the large background coming from down- 
going muons Q. Fig. |l] shows three different topologies of neutrino events: up throughgoing muons, 
semicontained upgoing muons (IU) and up stopping muons+semicontained downgoing muons (UGS + 
ID) @. The parent i/„ energy spectra for the three event topologies were computed with Monte Carlo 
methods. The number of events measured and expected for the three topologies is given in Table § 
Sources of background and systematic effects were studied in detail in || and found to be negligible. 
High energy data. The up throughgoing muons come from v^s interacting in the rock below the 
detector; the v^s have a median energy E v ~ 50 GeV. The muons with > 1 GeV cross the whole 
detector; the ToF information provided by the scintillation counters allows the determination of their 
direction (versus). 

Fig. ^ shows the zenith angle distribution of the measured flux of up throughgoing muons with energy 
^,1 GeV; the Monte Carlo expectation for no oscillations is shown as a dashed line, and for a — » v T 
oscillated flux with sin 2 29 = 1 and Am 2 = 0.0025 eV 2 , as a solid line. The data have been compared with 
Monte Carlo simulations, using the neutrino flux computed by the Bartol group and the cross sections 
for neutrino interactions calculated using the GRV94 parton distribution; the propagation of muons to 
the detector was done using the energy loss calculation by Lohmann et al. The total theoretical 

uncertainty on the expected muon flux, adding in quadrature the errors from neutrino flux, cross section 
and muon propagation, is 17 %; it is mainly a scale error that doesn't change the shape of the angular 
distribution. The ratio of the observed number of events to the expectation without oscillations is given 
in Table 0. 

The independent probabilities for obtaining the number of events and the angular distribution ob- 
served have been calculated for various parameter values. The value of Am 2 obtained from the shape 
of the angular distribution is equal to the value needed to obtain the observed reduction in the number 
of events; for Vn v r oscillations the maximum probability is 57 %, with best parameters Am 2 = 
0.0025 eV 2 , sin 2 29 = 1. The probability for no-oscillations is 0.4 %. The probability for za, — > v steri i e 
oscillations is 0.15%; combining these probabilities with the measured ratio of horizontal/ vertical muons 
(Fig. ^ we conclude that — > v ster n e oscillations are disfavoured compared to — > v T . Fig. ^ shows 
the 90 % CL regions for — *• u T , computed according to ref. Q. 
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Figure 2: The black points are the measured fluxes of the up throughgoing muons with > 1 GeV plotted vs. 
zenith angle 6. The dashed line is the expectation for no oscillations; it has a 17 % scale uncertainty. The solid 
line is the fit to an oscillated muon flux, obtaining maximum mixing and Am 2 = 0.0025 eV 2 . 
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Figure 3: The flux ratio vertical/horizontal for upthroughgoing muons. The measured point is plotted at a Ar 
around the minimum of y 2 for — > v T oscillations. 
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Figure 4: Confidence regions for i/ M —* v T oscillations at the 90 % CL calculated according to ||, obtained from 
MACRO high energy and low energy data. 
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Table 1: Event summary for the MACRO atmospheric neutrino flux analyses. The ratios R = Data/MC are 
relative to MC expectations assuming no oscillations (column 3). 



Low energy data. The Upgoing semicontained muons come from interactions inside the lower 
apparatus. Since two scintillation counters are intercepted, the ToF is applied to identify the upward 
going muons (Fig. |l|). The average parent neutrino energy for these events is 4.2 GeV. If the atmospheric 
neutrino anomalies are the results of — > v r oscillations with maximum mixing and Am 2 between 10~ 3 
and 10~ 2 eV 2 one expects a reduction of about a factor of two in the flux of these events, without any 
distortion in the shape of the angular distribution. This is what is observed in Fig. || 

The up stopping muons are due to external interactions yielding upgoing muon tracks stopping in 
the detector; the semicontained downgoing muons are due to induced downgoing tracks with vertices 
in the lower MACRO (Fig. [j]). The events were selected by means of topological criteria; the lack of time 
information prevents to distinguish the two sub samples, for which an almost equal number of events 
is expected; the average neutrino energy for these events is ~ 3.5 GeV. In case of oscillations with the 
quoted parameters, a reduction in the flux of the up stopping events as the semicontained upgoing muons 
is expected. No reduction is instead expected for the semicontained downgoing events (coming from 
neutrinos with path lengths of ~ 10 km). 

The MC simulation for the low energy data uses the Bartol neutrino flux and the neutrino low energy 
cross sections of ref. 0. The number of events and the angular distributions are compared with the 
predictions in Fig. ^| and Table [l]. The measured low energy data show a uniform deficit with respect 
to the predictions; there is good agreement with the predictions based on neutrino oscillations using the 
parameters obtained from up throughgoing muons. The allowed region in the Am 2 , sin 2 29 plane is shown 
in Fig. f| 
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Figure 5: Measured and expected number of low energy muon neutrino events versus zenith angle. Top graph: up 
stopping plus down semicontained; bottom graph: up semicontained. The solid lines are the predictions without 
oscillations; the dotted lines are the predictions assuming neutrino oscillations with the parameters obtained from 
the up throughgoing sample. 



Using the double ratio R — (Data/MC)iu/(Data/MC)m+uGS between data and MC of the two 
low energy data sets, the theoretical uncertainties on neutrino flux and cross sections almost disappear 
(a residual 5 % uncertainty remains due to the small differences between the energy spectra of the 
two samples). The average value of the double ratio over the measured zenith angle distribution is 
R ~ 0.59 ± 0.07stat. R = 0.75 is expected in case of no oscillations, R — 0.58 in case of — > v T 
oscillations. 



3 High Energy Muon Neutrino Astronomy 

The excellent angular resolution of the detector allows a sensitive search for up-going muons produced by 
neutrinos coming from celestial sources; in this case the atmospheric neutrinos are the main background. 
The pointing capability of MACRO was demonstrated by the observed " shadows" of the Moon and of 
the Sun, which produce a "shield" to the cosmic rays. We used a sample of 45 x 10 6 downgoing muons, 
looking at the bidimensional density of the events around the center of the Moon and of the Sun 0. 
The density of events shows a depletion with a statistical significance of 5.5 a for the Moon and of 4.2 a 
for the Sun. The slight displacements of the maximum deficits are consistent with the deflection of the 
primary protons due to the geomagnetic field for the Moon and to the combined effect of the magnetic 
field of the Sun and of the geomagnetic field. The angular resolution, corresponding to a cone angle 
containing 68 % of the events, is 0.9° ±0.3°; it is obtained from double muon events. An excess of events 
over the atmospheric v background was searched for around the positions of known sources in 3° (half 
width) angular bins. This value was chosen so as to take into account the angular smearing produced by 
the multiple scattering in the rock below the detector and by the energy-integrated angular distribution 
of the scattered muon, with respect to the neutrino direction. A total of 1197 events was used in this 
search, see Fig. |(| No excess was observed; the 90 % CL limits on the neutrino fluxes from specific 
celestial sources are at the level of ~ 10~ 6 cm -2 s -1 , see Table |§. 

We also searched for time coincidences of our upgoing muons with 7-ray bursts as given in the BATSE 
3B and 4B catalogues, for the period April 91 - October 99 M. No statistically significant time correlation 
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Table 2: High energy neutrino astronomy: muon and neutrino flux limits (90 % CL) for selected sources calculated 
according to the prescription in Q|. Previous best limits (B is for Baksan, I is for 1MB) were computed using the 
classical Poissonian method. 
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Figure 6: High energy neutrino astronomy. Upgoing muon event distribution in equatorial coordinates 
(1197 events). 



G 




m„ (GeV) 



Figure 7: Upward-going muon flux vs neutralino mass m x for Efj_ = 1 GeV from the Earth. Each dot 
is obtained varying model parameters. The solid line is the MACRO flux limit (90% CL); the line for 
the no-oscillation hypothesis is essentially indistinguishable in the log scale from the one for the — > v T 
oscillation hypothesis (it could be about two times lower). The open circles indicate models excluded by 
direct measurements (particularly the DAMA/Nal experiment) and assume a local dark matter density 
of 0.5 GeV cm- 3 . 



was found. 

4 Indirect Searches for WIMPs 

Weakly Interacting Massive Particles (WIMPs) could be part of the galactic dark matter; they could 
be intercepted by celestial bodies, slowed down and trapped in their centers. WIMPs and anti- WIMPs 
could annihilate and yield pions, kaons and then neutrinos of GeV or TeV energy coming from a small 
angular window from the celestial body centers. The 90 % CL limit for the flux from the Earth center is ~ 
0.8T0" 14 cm -2 s _1 for a 10° cone around the vertical. The limit from the Sun is ~ 1.4x 10~ 14 cm -2 s . 
If the WIMPs are identified with the lowest mass neutralinos, the MACRO limit constrains the stable 
neutralino mass, following the models of Bottino et al., see Fig. [?] S. 

5 Neutrinos from Stellar Gravitational Collapses 

A stellar gravitational collapse (GC) is expected to produce a large burst of all types of neutrinos and 
antineutrinos with energies of 7 — 30 MeV for a duration of ~ 10 s. The z7 e 's can be detected via the 
reaction v e + p — > n + e + in the liquid scintillator. About 100 -j- 150 v e events would be detected 
in our 580 t of scintillator for a stellar collapse at the center of our Galaxy. We use two electronic 
systems for detecting P e 's from stellar gravitational collapses ||. Both systems have an energy threshold 
of ~ 7 MeV and record pulse shape, charge and timing informations. Immediately after a trigger, the 
dedicated system lowers its threshold to about 1 MeV, for a duration of 800 /i s, in order to detect (with 
~ 25 % efficiency) the 2.2 MeV 7 released in the reaction n + p — > d + 72.2 MeV induced by the neutron 
produced in the primary process. A redundant supernova alarm system is in operation: we have defined 
a general procedure to alert the physics and astrophysics communities in case of an interesting alarm 
[||. A procedure to link the various supernovae observatories around the world (MACRO, SuperK, LVD, 



7 



SNO, etc.) was set up. Our live-time fraction in the last four years was ~ 97.5 %. No stellar gravitational 
collapses were observed in our Galaxy since 1989. 



6 Conclusions 

We presented new MACRO data on upgoing muons of low and high energy. These data are relevant for 
the study of neutrino oscillations, in particular — > v T ]To| | . 

We discussed the limits obtained from the searches for astrophysical sources of H.E. neutrinos, for 
WIMPs and for low energy neutrinos from gravitational stellar collapses. 
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